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Abstract 

In this paper I show that any mth-degree polynomial function of the elements of the density 
matrix p can be determined by finding the expectation value of an observable on m copies of p, 
without performing state tomography. Since a circuit exists which can approximate the measure- 
ment of any observable, in principle one can find a circuit which will estimate any such polynomial 
function by averaging over many runs. I construct some simple examples and compare these results 
to existing procedures. 
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I. INTRODUCTION 



Much of the practice of quantum information theory revolves around calculations of 
various information measures. These are functions of the quantum state p which quantify 
properties of interest, such as purity, entanglement, and distinguishability. An important 
subset of these are polynomial functions. Some polynomial quantities of potential interest 
are unitary invariants such as Kempe's invariant and the S-tangle , and measures of 

entanglement such as the Q measure of Meyer and Wallach 0, ^. Some other quantities, 
which are not polynomials themselves, can either be approximated by polynomials (such as 
the quadratic approximation of the von Neumann entropy) or are simple functions of such 
polynomials, such as the concurrence and the negativity of the partial transpose 0,0,0]. (I 
state what is meant by a polynomial function of p more precisely below.) 

A difficulty in estimating polynomial functions of the state is the linearity of quantum 
mechanics. Given only one copy of a system in a particular quantum state p, any mea- 
surement can only depend linearly on p. If it is possible to repeatedly prepare systems in 
the state p, we can overcome this obstacle by performing quantum state tomography P|: by 
estimating the expectation values in the state p of a set of observables, one can construct a 
classical description of the state; given that description, obviously any function of p can be 
calculated. However, many measurements are needed, of many different observables, since 
the number of parameters to be estimated grows like d"^ where d is the dimension of the 
system's Hilbert space. 

Recently there has been considerable interest in the problem of determining polynomial 
functions of quantum states (or simple functions of such polynomials); and a number o" 
ingenio us proc edures have been put forward for measuring particular functions 
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17, 



18|. These procedures try, by various clever techniques, to avoid 



estimating all the parameters needed for full state tomography; rather, they either estimate 
the polynomial function directly, or estimate a smaller set of parameters from which the 
polynomial can be calculated. They do this by performing measurements on several copies 
of p at once. 

In this paper I demonstrate that for any mth-degree polynomial / in the elements of the 
density matrix p (on a Hilbert space H), there are observables Of and O'j on Ti®"^ such that 
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the expectation value of 0/ + iO'j on m copies of p equals /: 

/ = TT{{df + .6})p«™} = {Of) + z{d'f). (1) 

II. POLYNOMIAL FUNCTIONS 

What do I mean by a polynomial function on the elements of p? Let us suppose that 
{l^)} is a basis for the Hilbert space Ti, where < i < d. Then we can write p in terms of 
that basis. 

d-l 
i,j=Q 

A polynomial function f{p) of degree m can be written 

f{p)= E ^^.n■..^^J^P^.nP^.n■■■P^r.Jr., (3) 

where the {ci^ji...imjm} arbitrary complex constants. There could also be terms of lower- 
order than m; but we can exploit the fact that Trp = pa = 1 to write them all as sums 
of terms of order m, so the form Q is general. This is easy to see: consider a term of order 
k < m, proportional to pj^j^ ■ ■ ■ Pi^.j^.- We can rewrite this as a sum of terms of order m by 
multiplying it by m — A; factors of Trp: 

ikjk / J Phjl ' ' ' PikikPik+lh+l 

■ -P (4) 

General polynomials of form (jS)) will not be invariant under a change of basis (though 
some polynomials of interest are invariant). However, since a basis change is a linear trans- 
formation, such a transformation takes mth-degree polynomials to mth-degree polynomials. 

Suppose we have m copies of the system in state p; these copies have the joint state 

P®"'= Aiii---Pw,nKl)Ol|®---® (5) 

Consider one term of the polynomial function /(p) given by 0, proportional to 
Phjl ■ ■ ■ pimjm- We can write down an operator 

Kji-imjm = ® ■ ■ ■ ® \jm){im\, (6) 

such that 

Tr{y4j^j^...i^j^p® } = Pi^j^ ■ ■ ' PimjmJ C^) 



as we see by substituting © and © in (|7j). Since the full polynomial /(p) is a linear 
combination of such terms, and the trace is a linear operation, we can find an operator Af 

~ ^ ] ^iljl---imjm^iljl--imjmy (S) 

h,ji,--;id,jd 

such that 

/(P) = (Af) = Tr{Ajp^-}. (9) 

While this operator obviously exists, it will in general not be an observable — that will 
only be the case if ^4/ is also Hermitian, Af = A^ However, we can trivially find a pair of 
operators which are observables. 

Of = (i/ + 4)/2, 

O'f = -z{Af - 4)/2, (10) 

such that {Af) = {Of) + i{d'f). 

These operators will not be unique, in general. If A/ is an operator as defined above, then 
the operator obtained by permuting the order of the m systems will also give the polynomial. 
So will any linear combination of such permutations, so long as the coefficients of the linear 
combination sum to 1. This freedom reflects the permutation symmetry of the state p®"* 



III. MEASURING POLYNOMIAL FUNCTIONS 



If we have the ability to measure arbitrary observables, then we are done; we repeatedly 
prepare p®™ and measure Of and O'f until we have a good estimate of their expectation 
values. If we do not have the ability to measure arbitrary observables, the fact that an 
observable exists whose expectation value yields the polynomial function is not very useful 
in itself. However, even if we cannot measure arbitrary observables, we can design a circuit 
which will simulate such a measurement to any desired degree of accuracy. 

The proof of this is trivial. Let the observable Of have a basis of eigenstates \(j)j) with 
corresponding eigenvalues oj (which are not necessarily distinct). Let \j) be a standard basis 
which we are capable of measuring. (For instance, in a system of q-bits this could be the 
computational basis.) Then we can deflne a unitary transformation 

Uf = J2\j)(^^\^ (11) 
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which rotates the observable 0/ to be diagonal in the standard basis: 

Of = Y.''M^)^^^^ UfOfU\ = Y.^,\3){3\. (12) 

j 3 

Clearly this new observable IfjOfUj is straightforward to measure: we just measure the 
state in the standard basis and give a value Oj to outcome j. This gives us an expression for 
the expectation value of Of, 

Tr{dfp^'-} = TT{{UfdfU}){Ufp^-U})} 

= E^^01(^/P^™f>))|j)- (13) 

j 

All we need then is to be able to find a circuit which approximates the unitary transformation 
Uf. It is well-known that a circuit exists to approximate any unitary transformation to any 
desired degree of accuracy Our procedure now is as follows: repeatedly prepare m 

copies of the state p, and carry out the circuit which approximates Uf. Measure in the 
standard basis, and weight each outcome j with the eigenvalue oj. Do the same for the 
observable O'f. We see, therefore, that in principle it is always possible to find a circuit 
which enables one to measure any polynomial function of p. 



IV. FINDING CIRCUITS 



This proof of principle is very far from being a proof that such a measurement is prac- 
tical. While such a circuit always exists, there is no guarantee that it will achieve the 
desired unitary transformation efficiently Indeed, most unitary transformations cannot be 
approximated efficiently j^]. A protocol of this nature would only be worthwhile if it were 
substantially more efficient than performing quantum state tomography on p and then cal- 
culating the function directly. This requires comparing the cost of doing more measurements 
(in the case of state tomography) to the cost of preparing a more complicated state (p®"^ 
instead of p) and performing some (possibly large) number of extra gates. 

Even if it is possible, one must in most cases diagonalize a pair of complicated observables 
Of and O'f and find the corresponding eigenvalues and eigenvectors. This may be a nontrivial 
task; if the system is of dimension d, then Of is a. d"^ x d"^ Hermitian matrix. However, 
for some applications this might be acceptable; the eigenvectors and eigenvalues could be 
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calculated "off-line," and then used to design a circuit which could be applied to many 
different states p. There is no known efficient algorithm for finding the simplest circuit that 
produces a given unitary transformation. Indeed, that is almost certain a computationally 
intractable problem in itself. But in some cases, it may be possible to find a simple circuit to 
do the job. Polynomials with a great deal of symmetry may simplify the circuits considerably. 

Let us consider a simple example. The first, and undoubtedly most immediately useful, 
is calculating the trace of 

This quantity gives a useful measure of the purity of a state, with Tr{p^} = 1 for a pure 
state and Tr{p^} = 1/d for the maximally mixed state; this is easier to calculate than 
the von Neumann entropy. For a bipartite pure state |\Ef), Tr{p^} gives a measure of the 
two subsystems' entanglement, where pa = Tr5{|^)(\l'|} is the reduced density matrix for 
subsystem A. From ()14|1 it is obvious that the correct observable is the swap operator: 

0/ = $^K)(j|®|j)(^|. (15) 

This is a special case of the widely known fact that 

Tr{p™} = Tr{5p®"^} (16) 

where S is the cyclic shift operator: Sl^pi) ® ■ ■ • ® \ipm) = I'ipm) ® li'i) ® • ■ ■ ® lipm-i). In the 
case m = 2, this is the pairwise swap, and is Hermitian as well as unitary, so that it is itself 
an observable. It has eigenvalues ±1, as any operator which is both Hermitian and unitary 
must; and if our system is a single q-bit, a good eigenbasis is the usual Bell basis: 

|0o) = (|00) + |ll))/v^, oo = l, 

101) = (|00)-|ll))/v^, 01 = 1, 

102) = (|01) + |10))/v^, 02 = 1, 

103) = (|01)-|10))/V2, 03 = -1. (17) 

Circuits for measuring the Bell states are widely known; Fig. 1 presents an example of 
such a circuit, consisting of a controlled-NOT gate, a Hadamard gate, and two single-bit 
meaurements in the computational basis. In many cases it would be much easier to measure 
such an observable than to do full state tomography even on a single bit. 
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FIG. 1: Circuit for measuring two qubits in the Bell basis. 
V. RELATIONSHIP TO OTHER APPROACHES 



Other circuits have been presented in the literature for calculating other polynomial 
functions of the state [3, ll^ 12, ll^ l^^- It is not difficult to see that these circuits give 
a method for estimating the expectation value of a suitable observable Of, as described in 
this paper. 

Let us look at this a bit more closely. Consider the circuit in Fig. 2, which is a version of 
the circuit presented in 16i]. In this circuit, a set of m systems is prepared in the state p*^™; 
an extra "control bit" is prepared in the state |0). This control bit undergoes a Hadamard 
gate; it then serves as the control for a controlled-^/ gate (or circuit) with the m systems, 
followed by another Hadamard. Finally, the observable Z = |0)(0| — |1)(1| is measured for 
the control bit. It is straightforward to show that the expectation value of this measurement 
is (Z) = (l/2)Tr{(?7 + f/t)p®'"}. If f/ = 5, for instance, then this would be equal to Tiip""}. 
The operator (1/2) (t/ + t/^) is Hermitian and therefore an observable. 

Why does this work? The values of the measurements can only be 1 and —1; for most 
observables these will not be equal to any eig envalue of (l/2)(f/ + W). In spite of this, 
the expectation values are exactly the same as if the observable (1/2) (f/ + f/"!") had been 
measured. To see what is going on, let us suppose for the moment that U has eigenstates 
\(f)j) with eigenvalues exp{i6j). We can always express the state of the m systems in terms 
of these eigenstates. 
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J2Rjj'\(pj){4>j'\- 



(18) 
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Suppose that we start the control bit and the m systems in the joint state = |O)|0j), 
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and carry out the circuit in Fig. 2. The resulting state will be 

m = eM^Ojm [cos(e,/2)|0) - ^sin(^^,/2)|l)] = exp(ze,/2)|7,)|0,). (19) 

The expectation {Z) in the control-bit state is cos^(6'j/2) — sin^(0j/2) = cos(^^j) which 
is the jth eigenvalue of the operator (l/2)(f/ + U'^). A similar result is obtained for each 
eigenstate This circuit is therefore equivalent to the unitary transformation 

E = ^e^^^/2exp(-i%X/2) ® (20) 
i 

So a different way of understanding this circuit sees the control bit as a target bit: if the m 
systems are in the jth eigenstate of U, the target bit is rotated about the x axis from |0) to 
a new state I'jj) whose Z expectation exactly matches the jth eigenvalue of the observable 
(l/2)(f/ + f/''"). If we now start the m systems in the state p®", the different eigenstates of 
U will contribute with their appropriate weights Rjj: 

Tr{(z®/ni(|o)(o|®pn^^} = J2^^rMZ\lj){<Pf\<Pj) 

if 

3 

= J]i?,,(0,|(l/2)(f/ + f/t)|0,) 
j 

= j;%'(0,'|(l/2)(f/ + f/t)|0,) 
jj' 

= Tr{(l/2)(f/ + f/^)p®™}. (21) 

Clearly, for any observable O whose eigenvalues oj lie in the range — 1 < < 1, we 
can find a unitary JJ such that O = (l/2)(f/ + [/^). For example, JJ = O + i^/ 1 — is 
such an operator. If the eigenvalues of O lie outside this range, we can construct a new 
observable 0/c where c > is larger than the absolute value of the largest eigenvalue 
of O. We can then find a unitary JJ such that (1/2) (f/ + JJ'^) = 0/c, and estimate the 
expectation value of 6/c by carrying out the circuit in Fig. 2. This is an alternative and in 
some cases much more efficient way of estimating a polynomial function of p; in particular, 
it only requires measurements of a sing le bit. Paz and Roncaglia Q have constructed a 
programmable circuit of this type to estimate the expectation value of any observable. This 
would obviously be sufficient to estimate any polynomial function of the state p. 
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FIG. 2: Circuit for measuring {\/2){U + W); after performing this circuit, the expectation value 
of Z for the first bit is equal to the expectation value of (l/2)(^7 + [/^) in the initial state of the 
other m systems. 



VI. CONCLUSIONS 



The ability to measure arbitrary m-system observables and determine their expectation 
values allows the estimation of arbitrary mth-degree polynomials of the state, without having 
to perform quantum state tomography. This is an enormous simplification, in principle; in 
some cases, it may give a considerable simplification in practice. Whatever the practicalities, 
this equivalence is a useful fact, which may make possible improved protocols for estimating 
entanglement measures, polynomial approximations to the entropy, or other functions of the 
state. 
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